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III.l AERODYNAMICS 


The aerodyneunic subroutines are concerned with forces and mcnnents as 
they vary with flight condition and attitude. The main progreun calls four 
programs which could broadly be termed aerodynamic in nature, as listed in 
Table I I I. 1.1. These programs are of lengths indicated in the parentheses 
and call two others, CLIFT and DRAG, and this volume of the documentation 
will be concerned with descriptions of all of these seven subroutines. 

III. 1.1 Subroutine AERO - 
Equivalent Flat Plate and Wetted Areas 

This subroutine is called by CTA.SR, and is principally concerned with 
a set of drag coefficients which are found by combining geometric descriptors 
with experimentally determined aerodynan\ics . 

Drag is assumed to be composed of profile drag, lift-induced drag, aiui 
compressibility drag. Profile drag is the most complex portion, and for each 
airplane component the corresfxinding drag coefficient takes the form 

ACpp = Cp F(Re)(S/S^) (II 1.1.1) 

where 

= geometry-dependent form factor 
K 

C_ = Mach number dependent skin friction coefficient at a reference 
F 

7 

Reynolds number of 10 

F(Re)“ Reynolds number correlation factor 

S/S * ratio oi cemponent area to wing reference area (fuselage and 
w 

nacelle wetted areas and lifting surface planform areas) . 



aero - COI-IPUTES COMPONENT EQUIVALENT FLAT PLATE AND WETTED AREAS AND PROFILE DRAG 
(130) 


AEROUT - PRINTS AND PLOTS LOW AND HIGH SPEED DRAG POLARS (170) 


CLIFT - DETERMINES LIFT COEFFICIENT OR ANGLE-OF-ATTACK (40) 


DRAG - DETERMINES DRAG COEFFICIENT (90) 


FLAPS - DETERMINES MAXIMUM LIFT COEFFICIENT AND DRAG INCREMENT FOR VARIOUS 


FLAP TYPES AND FLAP SETTINGS (350) 


CTAER - DETERMINES REQUIRED LIFT COEFFICIENT AND DRAG COEFFICIENT (50) 
IN CRUISE FLIGHT 


TABLE III. 1.1 - SUBROUTINE STRUCTURE OF AERODYNAMIC CALCULATIONS 




III. 1.1.1 Form Factors. — Computation of the form factors begins with 


the equation defining the average thickness to chord ratio for the exposed 
wing. 


TC 


(TC - S /B(TC - TC )1 * T S TC 

* R WF^ * R T ^ ERM LM T 

T + S ~ 

ERM LM 


(III. 1.2) 


where 


ERM 


1 - (1-S^. 


TC„,TC ■ thickness to chord ratios (root and tip) 
R T 

S ■ dicuneter of fuselage at wing, ft 

WF 

B « wing span, ft 

S ■ taper ratio (tip chord/root chord) 

LM 


In terms of this dimensionless parameter, the wing form factor is 


C^ - 1.03(2. + 4. TC + 240. TC ] 
ICW 


(III. 1.3) 


Other form factors may also be input and the computed value is overridden in 
GASP by a positive input value in this and subsequent form factor equations. 
The corresponding factor for the vertical tail has the same form. 


'KVT 


2. + 4. 


TC ^ + 240. TC,^ 
VT VT 


(III. 1.4) 


The factor for the horizontal tail is modified by its location (low or high) 
on the vertical tail; i.e., 

W “ * ♦ 4.TC^ ♦ 240.TC^1 (III. 1.5! 

where “ 0 for a low-tail geometry and ■ 1 for a T-tail geometry. 
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The form factors for the nacelles and tip tanks (if any) are respectively 


^KN " * .35/(ELj/DBAR^)J (III. 1.6) 

2Uld 

^KTP ■ ^ .0025(AXIS/BXIS) ♦ 60./(AXXS/BXXS) ^ (III. 1.7) 

where 

EL ■ length of nacelle, ft 
N 

DBAR^ " diameter of nacelle, ft 
AXIS « length of tip tank, ft 
BXIS - diameter of tip tank, ft 


For the fuselage, the corresponding factor is 


where 


^KF ’ *^‘^25(Z^/S^) + 60./(Zj^/S^^)^l 

(III. 1.8) 


2 •» fuselage length, ft. 

L*F 


III. 1.1. 2 Skin Fricwion Coefficient . — The skin friction coefficient is 

a function of both Mach number and Reynolds number (Figure III. 1.1). It is 

7 

initially computed at a reference Reynolds number of 10 . At this Reynolds 
number its Mach number dependence is approximated by: 


C 


FIN 


.002944 

.00299 

.0027; 


.000176 E..; E.. < .4 

M M 

.0003125 E„; .4 .< E. < .92 


(III. 1.9) 


where E,, is the Mach number. This function varies nearly continuously over 
M 

the inter/al .0027 to .002944. 
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T'lII 

SKIN FRICTION COEFFICIENT 


FIGURE III. 1.1 REFERENCE SKIN FRICTION COEFFICIENT 





Conponent Mach-dependent skin-friction coefficients at the reference 
Reynolds number are next defined as equal to this definition is effec- 

tive for wing, nacelle, fuF^^lage, verticAl tail, horizontal tail, and tip tank, 
and they are respectively denoted by and 

^DTIP' 


III.1.1.3 Reynolds Number Correction . — The skin friction coefficient of 
Section III. 1.2 is determined at a reference Reynolds number of 10^. Effects 
of ccxnponent Reynolds number of the skin friction coefficient axe incorporated 
through the Prandtl-Schlichting turbulent flat plate skin friction equation. 
This results in the Reynolds number correction factor taking the general form 
of 


' 4‘:S!!E2^ . 


- 2.6 


(III. 1.10) 


'^Re“10 

The component Reynolds number is determined from the Reynolds number per foot 
(RELI) and a characteristic length. RELI is calculated at the design cruise 
conditions (HNCRU, EMCRU) and the characteristic length is usually taken in 
the direction of the airflow. 

In the case of the fuselage, we have 


ELp - fuselage length 

ELpcc “ fuselage length excluding tail boom 
and thus the fuselage Reynolds number is 


RF. 


“li Wo’'! 


*^I ^^FFC' *^C0NFg“^ 


(III. 1.11) 
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Similarly, other Reynolds numbers are fotutd as 


RE^„ ■ RE. _ • EL_„ (fuselage boom) 

BM LI BM 

Hv • RE. _ • EL., (nacelle) 

LX N 

RE^ •» RE^j • CBARy,j, (vertical tail) 

REj^ ■ RE^ j • (horizontal tail) 

*^TP " "^Ll * 


(111. 1.12) 

(111. 1.13) 

(111.1.14) 

(111.1.15) 

(111. 1.16) 


where EL...* .... AXIS, are characteristic lengths in the streamwise direction 
BM 

of the indicated ai.'craft components. 


The Reynolds numbr.:. correction factor for each component then becomes 


(wing) 

FFj^ - ^^10 (fuselage) (III. 1.17) 

“ ^^^10 ^ ^ (fuselage boom, if (III. 1.18) 

FNj^ - ILOG^q(RE^)/ 7]"^*^ (nacelle) (III. 1.19) 

FVTj^ - [L0 Gj^q(RE^)/ 71“^*^ (vertical tail) (III. 1.20) 

- lLOG_(RF. )/7l"^'^ (horizontal tail) (III. 1.21) 

KL XU HI 

FLTj^ - [LOGj^Q(RE^p)/7r^*^ (tip tanJts) (III. 1.22) 

III. 1.1. 4 Profile Drag Coefficient . — The profile drag coefficients of 
all aircraft components except the wing are obtained by summation. 


III.l 
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“s - «*p ®ri "re 

* CKyr ^VTI '^RE ®VT^®W 

* “ht “hti ^“"re 

* <\ “»I ™RE 


+ CK C FLT S /S ♦ DEL 

DTIP RE tip' W CD 

+ CK„„ S^«/S„ ♦ GR^„ (III. 1.23) 

BM DTIP RE WBn W CD 

where DEL is a fuselage incremental F which may be input; DEL„ is an incremental 
FE e CD 

input drac coefficient, and is the exposed landing gear drag coefficient. 

This last term i- included here only for fixed landing gear configurations 

^ O);for retractable gear configurations (I-„.„ • 0) it is added in 
GEAR GEAR 

subroutine DRAG. 

Landing gear drag can be input as an equivalent flat plate area (GRl-'E) 
or the default value is ccmiputed as a function of gross weight, which is 


GR 


CD 


(.0032/S^)W^ 


.8 


(III. 1.24) 


Wing profile drag is handled in subroutine DRAG; however, the product of two 
of the factors used in DRAG to compute the wing profile drag axe computed in 
AERO; 

SA^ - CKW * FW„ (III. 1.25) 
6 RE 
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III. 1.1. 5 Wetted Areas aM Flat Plate Areas . — The next set of compu- 
tations in the subroutine deal with ccnnponent wetted areas euid equivalent flat 
plate areas. The wetted area expressions are all straightforward except for 
the wing which is 

®wf 

SWET„ - 2S„ - ^ yB" (III. 1.26) 

where S * wing area, sq ft 
w 

^WF * fuselage diameter at wing, ft 
S. „ = taper ratio of wing 

B « wing spcui, ft 

The wing wetted area of Equation (I II. 1.26) is derived in Appendix II lA. 

The six component flat plate areas are listed next; for example, for the 
wing (W) , 

% ■ “wi “w (III. 1.27) 

where the four terms on the right have been defined. This function has the 
units of area and is proportional to the drag. The remaining five components 
are those associated with fuselage, vertical tail, horizontal tail, nacelles 
and tip tanks. The sum of these coefficients is tnen expressed as 

FE * FE„ + FE„ + FE,^ + FE„„ + FE^, + FE^„ + DLTAFE (III. 1.28) 

W F VT HT N TP 

where the last term is 

DLTAFE - DEL^^ * S„ (III, 1.29) 

CD W 
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III. 1.1. 6 Concluding Computations . — The Oswald efficiency factor, used 


to compute induced drag, is computed by 

SEE - [1.035 + 1.19 AR(CDpyCOS^(RLM^^) + (III. 1.30) 

where CD_^ • (FE-FE )/S , wing-free profile drag coefficient 
PO w w 

CDp^ - FE^S^ , wing profile drag coefficient 

2 

The coefficient multiplying C in the induced drag equation is next 

Ij 

defined as 

SA^ - l./(PI * SEE * AR) (III. 1.31) 

This coefficient is reduced if the aircraft carries tip tanks L use of the 
effective increase in wing aspect ratio; i.e., in ' ' xs case, 

SA^ = SA^/(l+.5 B^jg/B) (III. 1.32) 

where B has been defined in Equation (II I. 1.7) 

XIS 

For fixed landing gear (J^ FE is incremented by GRFE as given in 

GEAR 

Equations (III. 1.28) 

FE - FE + GRFE (III. 1.33) 

III. 1.1. 7 Co mpressibility Drag Parameters . — The wing geometry enters 
into the next calculations with the specification of the tangent of the quarter 
chord sweep angle, 

T = *^^^^*^2104) (III. 1.34) 
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Following this is the paraneter depending on the taper ratio S , 

LN 




(III. 1.35) 


and sweep «mgles 


DLMpS - ATAN(T - VAR(X^pg - .25)YAj^Qgl * RTOD 


(III. 1.36) 


(III. 1.37) 

where 

y ^ wing minimum pressure point in fraction of wing chord 


y m wing maximum thickness point in fraction of wing chord 

CTCMX 

. .35 for conventional airfoils 


.40 for supercritical airfoils 


These points enter into the computation of one of the empirical compressibility 
parameters. 

SA^ - -.33(.65 - X^pg)[l + . 0033(4. *DLMpg-3.*DLM^^)J (in.1.38) 
The sweep angle of the wing leading edge is then 


RLM, ^ - ATAN(T + YA_^_/AR] 
laE) IjEOD 


(III. 1.39) 


and this is follovied by the definition 


''k * >1*«LE05''“ * ^ 


(111.1.40) 


Three more of the compressibility drag parameters are then defined as follows: 


III.l 
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SAj^ - tl + .033(4.*DLMpg - 3.*DLM^jj)J 

* 11. - 1.4*TC - .06(1. - S^pg)J - .0368 (III. 1.41) 

SA3 - (1.5 - 2.*F^ > , J^^(RLM^)1 (III. 1.42) 

2Uld 

SA^ - .75 • TC (III. 1.43) 

III. 1.2 Subroutine AEROUT - 
Drag Polar Computation and Plots 

While this is an apparently long subroutine, the major part of it is 
concerneid with the graphical presentation of routine lift-drag relationships. 
Thus, low speed and cruise polars are computed by calls to the subroutines 
CLIFT and DRAG using as input the Mach number and certain gross geometric 
characteristics of the aircraft. The coefficients of lift and drag are then 
plotted with angle of attack as the independent parameter in the usual drag 
polar format. That is, either the lift coefficient is specified and the angle 
of attack and drag coefficient are found, or the angle of attack is specified, 
and both lift and drag coefficients are founds In either case, C is then 

Li 

plotted on the vertical axis as the dependent variable with C^ as the independent 
variable. 

If the input flag NPC is greater than 1, the aircraft is in a cruise 
configuration, and the computations between statements 1 and 50 deal with the 
determination of the drag coefficient as both lift coefficient and Mach number 
are independently varied. The outer loop varies C from .1 to .5, and at each 

1j 
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value the inner loop varies the Mach number from .5 to .85. For each 

combination (C , the subroutines CLIFT, DRAG return the quemtities 

L 

ACD(I, J) - drag coefficient (C^ , X„„) 

L EM 

ZLOD « lift to drag ratio 

CL^t.ph * lift curve slope, per rad. 

ALPH* ■ angle of attack, deg. 

If the input flag NPC * 2, the angle of attack is the independent 
variable at ten values from -2 degrees to 16 degrees, and Mach is constant at 
.15. For each value, the subroutines CLIFT and DRAG return lift and drag 
coefficients in the three flight conditions, 

I “1: low speed cruise configuration 

Ij 

I * 2: gear up take-off configuration 

L 

I =3: gear up approach configuration 

L 

The lift and drag cirrays are tlien 

ACLKI . K) - lift coefficient (I , ALPHA) 

L L 

ACDld , K) = drag coefficient (I , ALPHA) 

li li 

The final portion of the subroutine is concerned with the plotting of the 

drag polars. First, the drag coefficient is plotted versus Macii number for 

eight constant values of lift coefficient. These correspond to the criuse 

configuration. Finally, the lift drag polar is drawn at low Mach number in 

the three flight conditions defined by the indicator I . 

L 
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III. 1.2.1 Subroutine CLIFT - Wing-Alone Lift Coefficient . — This 
subrout le relates the lift coefficient and angle of attack of the wing alone. 
Either the lift coefficient or the angle of attack c^m be input, and the other 


is determined, for the appropriate cruise, take-off or approach, and landing 


configuration. The input flag takes one of six values, reflecting 

LIFT 

these various possibilities; i.e., 


I.,.,™ “ Is Compute C, at cruise from ALPHA 
LIFT L 

2: Compute C at take-off from ALPHA 

L 

» 3; Compute C at approach from ALPHA 

L 

= 4; Compute ALPHA at cruise from C 

L 

= 5: Compute ALPHA at take-off from C. 

L 

» 6: Compute ALPHA at approach from C. 

L 

The lift coefficient and angle of attack are linearly related by the 

slope CL . Lift increments are those due to flap deflection (DEL ) 

^ ALPH CL 

lud to a multiplier for ground-effect (KCL ) . These are set equal to 0 and 

GE 

1, respectively, at the start of the subroutine. The first computation is 
that for the lift-curve slope, 

•• C^C^ (III. 1.44) 

where 

“ AR/12 COS(RLM^^))^ (III. 1.45) 

C^ * 1 - [EM * COS(RLM^^)]^ (III. 1.46) 

AR =■ aspect ratio of wing 
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RLM^^ ■ quarter chord sweep angle, rad 
EM ■ flight Mach number 

In the cruise configiuration, neither flap effect nor ground effect is 

significant, and I, ___ ■ 1 or 4. When I. " 1, ALPHA is iiput, and the 

LIrT LIFT 

cruise lift coefficient is given by 

- CL^^(ALPHA - ALPHLO) (. 0174533) + (III. 1.47) 

where 

ALPHA ■■ angle of attac)c, deg (input) 

ALPHLO ■ zero lift angle of attac)t, deg (input) 

This will be the lift coefficient returned unless CL^^ exceeds CLj^^ = 2. 

In this case, the maximum angle of attack is computed as 

ALPHA = ALPHLO + (CL^^ - DEL^^) /(. 01745 ♦ CL^^^) (IH.1.4H) 

and the lift coefficient returned is 

S " 01745) (ALPHA - ALPHLO) + DEL^^ (III. 1.49) 

1, Tc... ~ cruise lift coefficient is input, the program 

LIFT 

cOTtputes, 

^■■OGE ■ * “‘■CL (III. 1.50) 

where the right side reduces to CL, since KCL__ * 1 and DEL„ “ 0. 

G£ CL 

Then# if this is less t!ian CL the corresponding angle of attack is given 

MAX 

by 

III.l 
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ALPHA - ALPHLO ♦ (CL - DEL. ) / (1CCL__ CL.,_„ .01745) (III. 1.51) 

CL GE aLPH 

If “ 2 or 5, the tiOce-off configuration requires consideration of 

both flaps ^md ground effect 2 Uid of the parameter HOB. This is the altitude 
of the aircraft mac as measured in wing spans « and it enters into the compu- 
tation of the factor 


KCL^g “ 1.005 + C^(. 00211 - .0003(AR - 3)) (III. 1.52) 

where 

Cj - exp (5.2(1 - HOB)l (III. 1.53) 

For typical aspect ratios, this equation shows that CKL__ vuies between 

GE 

1.005 and 1.5, approximately, with the larger values occurring at low alti- 
tude (HOB = 0) . 

If IfTr-r ® angle of attack is input, the lift coefficient is 

liXf X 

again found using Equation (III. 1.47). otherwise *5, and the angle of 

LIFT 

attack is again computed according to Equation (III. 1.51), where DEL =DCL 

CL TO 

and ~ appropriate to the take-off configuration. 

If 3 or 6, as for the approach configuration, the computations 

use different values for the quemtities. 


DEL = DCL ^ 
CL LD 


CL = CLMX 
MAX LD 


but otherwise the same as for the take-off configuration. 


(111. 1.54) 

(111. 1.55) 
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III. 1.2. 2 Subroutine DRAG - Total Aircraft Drag Coefficient .— This 
subroutine is concerned with the total aircraft drag coefficient at a given 
flight condition and configuration. Profile drag, induced drag, and compressi- 
bility drag contributions are found, and ground effects, flap defleccion effects 
and landing geau: effects are included when appropriate. These flight conditions 
are distinguished by the calling parameter defined below. 

Plight Condition 

Cruise, gear up 
Take-off, gear down 
Landing, gear down 
Cruise, gear down 
Take-off flaps down, ge^u: up 
Landing flaps down, gear up 

7 Cruise, gear down 

8 Approach flaps, gear up 

Other input call parameters are the Mach number EM, dimensionless height above 
ground, HOB, and lift coefficient, C , and return quantities are drag divergence 

ia 

Mach number EM^ and drag coefficient C„. 

D D 

In the "clean" configuration (IDRAG = 1) , the aircraft drag coefficieiit 
is expressed as 

- SA_ + SA- CD,._ + CIC. DEL^^„ + SA_cf 
D 5 6 WI w CDM 7 L 

Profile Compressibility Induced 


^DRAG 

1 

2 

3 

4 

5 
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where 


SAg ■ profile drag of all components except the wing (excludes 
landing gear if retractable) (from AERO) 

* ™RE 

« wing form factor x Reynolds number correction factor 
(see AER>^) 

CD.., « wing s)<in friction coefficient 

W1 

DEL « compressibility term (see below) 

CDM 

SA^ ■ l/?TARe (from AERO) 

The wing profile drag may be expressed as a function of lift coefficient by 

multiplying the wing reference skin friction coefficient (CD ■ CFIN, from 

W X 

AERO) by an additional correction factor, COCDR, which is a function of lift 
coefficient. Values for CDCDR, along with corresponding table values of lift 
coefficient, CLS, must be input by the user. If CDCDR is not input, it is 
assumed to be 1.0 for all lift coefficients. 

The correction factor CDCDR represents the ratio of the actual clean wing 
F^rolile drag at some lift coefficient to the profile drag of the wing at the 
angle of attack for minimum profile drag. Typical values for CDCDR are illustra- 
ted in Figure II I. 1.2. Also illustrated in a set of values for CDCDR which would 
be representative of a laminar flow wing. Values less than 1.0 are appropriate 
in the "drag bucket" since the skin friction coefficient is computed assuming 
fully turbulent flow (see AERO) . 

The compressibility drag term is given by 
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10 (EM-EMjj)^, if EM EMp 

0* otherwise (II I. 1.56) 

where EM^ - SA^ ♦ SA^C^ + .00 (III. 1.57) 

•> drag divergence MacH nuniber 

The factors SA^ and SA^ are ccwiputed in AERO, ^ factor which may be 

input (as to shift the drag divergence Mach number of the airfoil. 

Figure (III. 1.3) 

For flaps and landing gear are down, and the drag coefficient 

is expressed as the sum of five terms, 



where 


and 


• SA_ + SA,CD„^ + del + DEL^„ 
D 5 6 WI CDF CDG 


+ SA_K - SIGMA * DEL^ ) /VDEL^ (111.1,58) 

' CDXGE L CL o 


DEL^^p ■ parasite drag coefficient increment due to trailing edge 
flaps, from FLAPS 

DELcjjg “ parasite drag increment due to landing gear, equal to 
as found in AERO 

SIGMA ■ induced drag correction factor, from FLAPS 

DEL - incremental lift due to flaps, from FLAPS 
CL 

VDEL * wing efficiency correction factor due to change in span- 
6 

wise lift caused by flap deflection, frcxn FLAPS 

2 1/2 

* 1.111 + 5.55 * HOB - [29.8(HOB + .02) + .817) (III. 1.59) 

COIGE 
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Ill 


SUPERCRITICAL 

ASSUMPTIONS 

(1) -AM^ = 0.08SCFAC 

(2) = 10 (M 
M 

(SAME FOR COt>n/ENTIONAL 
WITH LOWER Mj^) 



MACH NUMBER 

FIGURE 1 1 1. 1.3 - SCHEMATIC M. .1 EL OF SUPERCRITICAL AERODYNAMICS 



This last expression is to be used only ^hen in ground effect (HCB < .88, 
approximately) since the next line of the subroutine sets ^ 1, if this 
function is greater than 1. 

The final lines of DRAG involve the contribution of the nacelle drag 


c ■ c ♦ c 

D D DNAC 

where is defined in subroutine ENGS2 when KNAC ■ 1. 


(Ilia. 00) 


Ilia. 3 Subroutine CTAER ‘ Cruise Aerodynamics 

The cruise-drag computation, as indicated in Taole III. 1.1, is performed 
in the 50-card subroutine CTAER, which itself calls on AL.-l, CLIFT, and DRAG. 
Each of these subroutines has already been described in previous pages of 
this section, while CTAER itself is concerned only with details of the canpu- 
tation. 

The cruise drag is found as a function of cruise Mach number and altitude 
which are specified. The subroutine finds the angle of attack required to 
balance lift and weight. The drag (or thrust required) is found at this equi- 
librium flight condition. 

The cruise lift coefficient is given as a function of angle of attack in 
subroutine CLIFT, while the drag coefficient is obtained by calling aERO and 
DRAG, which combine several drag components as described in prior sections of 
this chapter. 

The cruise weight option is ''ecified by the indicator IDC, such that is 
IDC •= 1, the weight W is input through the UNIV common block. Otherwise, tin; 
ratio of cruise weight to design gross weight (RW^^^^^) and gross weight 
together imply W: 

W »= WG * (III.l.Gl) 

TO CRTO 
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At Statement 6, the subroutine TPALT provides veurious physical charac- 
teristics of the atmosphere as Sanctions of the gemetric altitude, 

N 

The static pressure, P , speed of sound, ASONO, and velocity V , are 

o o 

found in the required units, for the input values of altitude and Mach number, 
and the lift coefficient required for cruise is calculated as 




W/(.7 


®REF 


P 

o 


S 


2 

MN 


) 


III. 1.62) 


where S._ ic the Mv»' h number. The subroutine CLIFT then returns the angle 
MN 

of attack required in this cruise configuration. 

The subroutines AERO amd DRAG then return the corresponding drag coef- 
ficient, following which the cruiso lift and drag are found as 

' '^L (in. 1.63) 

DRG » ZL/ZLQV) (III. 1.64) 


where their ratio is ZLQD = C^/C„. 

L D 

I I I. 1.4 Subroutine FLAPS - 

Fla^i Aerodynamic, Geometric, and Weight Characteristics 
This subroutine is called by the main program and by two others (APPFLP 
and RGBAL) and it calls only the utility programs related to the current flight 
conditions and to interpolation of tabular data. The major portion of the 
subroutine is developed from a technical paper published by K. L. Sanders in 
May 1969, in the Journal of the Society of Aeronautical Weights Engineers. 

The purpose of the subroutine is the computation of aerodynamic, geometric, and 
weight characteristics of wing high-lift devices. A major output is the maximum 
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lift coefficient of the aircraft in the take-off or lemding configuration. 

Parasite drag incrcnicnts due to flap deflection ^u:e also founds as are flap 
areas and weights. 

Aerodynamic lift and drag data for a reference or baseline wing are known 
for various types of high-lift devices, and the program is chiefly concerned 
with tlie computation of correction factors, as applied to the reference wing 
configuration. Seven different flap types may be simulated as indicated by 
the input variable JFLTYP: 

JFLTYP « 1, plain 

“ 2, split 

* 3, single slotted 

■ 4, double slotted 

■ 5, triple slotted 

■ 6, Fowier 

* 7. double slotted Fowler 

The first third of the progreim listing is devoted to common card state- 
ments, comments and input data tabulations. This is followed by the GO TO 
statement depending on the input flag NPC. These options are briefly described 
as 

NPC ■ 2: Compute sensitivity factors leading to lift, 

drag characteristics for specific choice of 
flap 

NPC ■ 4: Compute flap lift and drag as flap deflection 

varied by subroutine APPFLF 
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III. 1.4.1 Ccxnputation of Airplane Geometry Lift and Drag Correction 


Factors . The basic standard wing has both flaps and leading edge slats, neither 
of which is deflected in the reference configuration. The rcferonce wing is 
assumed to have a thickness ratio of .10, and an aspect ratio of 12. It is 
unswept and untapered, and has a trailing edge flap chord ratio of .3 and a 
leading edge slat chord ratio of .15, both of which extend over the whole span. 
The reference Reynolds number and Mach numljer are 6 x 10 and 0 respectively. 

The maximum airplane lift coefficient is expressed as the sum of four terms. 


C • ^ R V V "f DCL V 

LMAX CLAMX L.V41 LAM2 MTE LAM3 ' 


. V 


LAMS 


^^MLE '^LAM9 


V IV V + DFL 

LAM12* LAM13 LAM14 CLF 


UII.1.65) 


where R “ input reference mciximum lift coefficient for basic wing 

CLMAX 

DCL “ input trailing edge flap lift increment for flap type at 
optimum flap deflection 

DCL “ input leading edge slat lift increment for leading edge typo 

at optimun> leading edge deflection 

DEL ■ fuselage lift increment 

CLF 


V , _ = sensitivity to variations from reference values for flaps 

LAMI 

and slats, I = 1 to 14 (correction factors) 


The first two correction factors relate to the wing aspect ratio and tldck- 

ness, and the next six (V, to V, correct for flap geometry and deflection 

LAM3 LAMO 

and for wing sweep. Similarly, V to V _ are used to account for 

LtAMX 4b 

variations in slat geometry and deflection, and the last two, and 

correct for Reynolds number and Mach number deviations from theii reference 
values . 

These correction factors are determined by successive calls to the inter- 
polation routine ITRLN and BIV. For example, for 
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CALL ITRLN(AAra«, ALAM1« AR, VXAMl, NARH) 


where AARW > input array of wing aspect ratios from 0 to 20, NARW in number 
ALAMl ■ corresponding input array of sensitivities of maximum lift 
coefficient to aspect ratio, NARW in number 
AR ■ actual input value of wing aspect ratio (COMMON/VNIV/) 

VLAMl » sensitivity of clean wing maximum lift coefficient to wing 
aspect ratio, AR 

NARW ’ number of data points in table AARW vx. ALAMl 
Subsequent calls to ITRLN and BT.V generate corrections for the following 
effects 


V 

LAM2 

V 

LAM3 

V 

LAM4 

V 

LAMS 

V 

LAM6 

V 

LAM7 

V 

LAMB 

V 

LAM9 


= sensitivity of clean winq C, to wing thickness to chord 

LMAX 


ratio 

= sensitivity of flap clean wing C 


LMAX 


to wing aspect ratio 


= sensitivity of flap clean wing C slope to wing thickness 
to chord ratio. This is a function of flap type 
= sensitivity of flap clean wing to wing flap to chord 

ratio (function of flap type) 

= sensitivity of flap clean wing C 
(function of flap type) 

“ sensitivity of flap clean wing 


LMAX 


to wing flap deflection 


to wing flap span 


sensitivity of flap clean wing C, to wing sweep angle 

LMAX 


o cos (SWPQC) 

» sensitivity of slat clean wing C 
= 6.65 • 

LEOC 


LMAX 


to slat chord 


(III. 1.66) 


(III. 1.67) 
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'^LAMlO 

V 

LAMll 

V 

LAM12 


- sensitivity of deem wing to slat deflection angle 

LNAX 

■ sensitivity of slat clean wing C. to slat span 

lmax 

“ sensitivity of slat clean wing to leading edge sweopback 

LMAX 

- cos^(SWPLE - 5/57.29) (III. 1.68) 


c is computed at the unaccelerated stall speed of the airplane corres- 
laNAX 

ponding to C . This is accomplished by an iteration on stall speed: at 

XiNAX 

each stall speed the Reynolds ntunber (V,...,.) and Mach number (V, correction 

LAM13 LAM14) 

factors are determined. Using these factors C is computed and based on 

XaMAX 

this C, a new stall speed is calculated. The iteration continues until the 
LMAX 

stall speed converges. 

The incremental lift due to trailing edge devices is given by: 


DCL = DCLMTE V, 

LAM3 


.V V V 

LAMS LAM13 LAM14 


Further computations deal with the modifications to the drag coefficient, as 
implied by the following factors: 


Vdeli “ sensitivity of flap to flap chord ratio (functions of 

flap size 

- = sensitivity of flap to flap angle, function of flap type 

DEL2 DniN 

Vdel 3 ” sensitivity of flap partial flap spam 

^ = sensitivity of flap to flap hinge line sweepbac)c 

DEL4 DrilN 

= cos(SWPp^^) (III. 1.69) 

c " scnsltivity of to fuselage width to span ratio, 

UEL^ DNXN 

= 1. - FWOB (III. 1.70) 

'^r^r.Tc * sensitivity of span efficiency factor to flap angle 
DELd 
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Sensitivities not followed by equations are determined by Interpolation of 
tabular input data, as was done for the lift coefficient incronents. The 
result of ct^nbining these expressions is the incremental drag coefficient 

DCD * ^*^OTE '^DELl * ' * '^DELS' parasitic drag 

coefficient, due tc trailing edge devices 

The empirical flap efficiency factors, SIGM^^ and SIGM^ are determined as 
functions of flap deflection for the take-off (IFLAP ■ 1) and landing (IFLAP 
» 2) configurations respectively. This relationship is illustrated in Figure 
IIT.1.4) 


III. 1.4. 2 Flap Weight Estimation . The flap weight is assumed to be a 
function of the flap type, flap area, and design flap airspeed. The flap 
area (ratioed to th<.) wing area) is computed irom the wing geometry: 


Slaps = “"''oc * ^ ‘ U - il - 1 


- FW_ - B^.,^^„]/(2 * S) 
OB OB ENGOB 


(III. 1.71) 


where Z 


=• root chord at wing centerline, ft. 


CF = ratio, flap chord to wing chord 

oc 

BTE^g = ratio, flap span to wing span 

= ratio, body dicuneter at quarter chord to wing span 
= ratio, total engine diameter to wing span, for wing mounted 

LNGOB 

engines 

The design airspeed of the flaps is assumed to be 


V “18V 
FLAP ■ STALL 


(III. 1.72) 


whore Vgr Ri.T. stall speed in knots corresponding to in the landing 


configuration. 
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Flap weights are represented as one of the following functions, depending 


on the flap type, Jp^TYP " ^ 


•^ELTYP “ 


”flAP “ (III. 1.73) 


*^FLTYF " 


2733 

.369 WC 2^, V* 

FLAP FLAP FLAP 


160 

FLAP — 


W, 


FLAP 


<“<=rLAP ''flap > ‘“ 


(III. 1.74) 


J * 3,4,5; Single, double or triple slotted flaps 

FLTYP 


'^FLAP* '^SlAP ®FLAP ^'^FLAP'^^°°'^ 


(III. 1.75) 


JpLTYp = 6, 7: Fowler or double slotted Fowler flaps 

si; (V^, ^„/100 . ) ^ FLAP 
FLAP FLAP FLAP FLAP N 


(III. 1.76) 


where WC„, = (.62, 1.0, .733, 1.2, 1.32, .633, .6781] = variation with 
FLAP 

Jfltyp flap weight factor or may be input 

S„, * S„, * S = flap area 
FLAP FLAPS ^ 

FLAP^, ■ 1. = number of flap segments per wing panel 
N 


The weight of the leading edge devices is given by the function: 

= 3.28 S^’^^ 

LED LE 


(III. 1.77) 
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where S is the leading edge slat area and the total weight of "high-lift 

laE 

devices" is the sum 

”hLDEV “ ”fLAP * ”lED (III. 1.78) 

111. 1.4. 3 Allowable Flap Reflection on Approach . The final section of FLAPS 

is used only diuing engine sizing when FLAPS is called by APPFLP with NPC = 4. 

FLAPS computes the incremental flap lift and drag coefficients and C. 

LMAX 

corresponding to an approach flap angle specified by APPFLP. Since engine 
sizing occurs only after the flap characteristics have already been determined, 
it is not necessary to recompute most of the correction factors (VLAM, VDEL) , 
and thus these sections of the program are bypassed when NPC ■ 4. Only the 
factors which are explicit functions of flap deflection are recomputed 
(VLAM6, VDEL2, VDEL6) . In addition, an empirical flap efficiency factor, 
SIGMAP, is determined from the flap angle specified by APPFLP. 

1 1 1. 1.4. 4 Subroutine APPFLP . This subroutine determines an approach flap 
deflection by calling FLAPS with NPC * 4 such that the corresponding approach 
lift coefficient is a specified fraction of the landing lift coefficient. The 
FAA regulation 25.121 requires that the approach stall speed cannot exceed 

1.1 times the landing stall speed, and this implies a specific flap deflection; 
i.e. , 

^LMXAP “ ^LMXLD'^^'^ * ^LMXLD 

where AP and LD denote "approach" and "landing" values, respectively. The 

landing value of C, has been previously computed by the program, and 
LMAX 

subroutine APPFLP returns the flap deflection DFLAP, as found by numerical 
intcrix)lation, that provides the appropriate approach C . 
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APPENDIX I 11. A 


WING WETTED AREA DERIVATION 



The exposed wing wetted area is 


SWET = 2S - 2S 

W W BLOCKED 


where it is assumed that the fuselage completely blocks the area shaded. Now 


BLOCKED 


/cr Cp \ 
CED “ I 2 j 


WF 


and for a wing with single taper 


Sup 

c = C X - C ) 

F R B R T 


substituting 


S =2 

BLOCKED 


[Cr+Cr-£^(Cr-Gj.) J 

B 




Now 


5 * * - lix ’ 


III.l 
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so that 


S ■2 

BLOCKED 


, c s 

3 WF 




a ♦ ^ (X-i) 


f c V 'i'*’ 


SwF 


c <> wr • - ~inrr .. % \ \ 

“blocked “ X+1 +X(X+1> -X ^ 2B ^ 


- 2 


^C S 
2 ^ WF 


il_Ly_ (1 + fwL ( 

■’ ^ 2B ‘ 


XH+X‘4X-\ 


^ X%X+1 2B 


therefore 


WET, 


W 


2S - 3CS - ^ 

^ X^X+1 


, WF , , 
1 + — ( ; 

2b 


-i)j 
- 1)J 

^ - 1 ) 
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The aerodynamics subroutines are called in the following order, in the 
typical application of program MAIN: 

1. FLAPS - Lift and drag dependence on wing geometry 

2. CTAER - Cruise lift and drag computation 

3. AEROUT - Drag polar variation with flight condition 

As described in earlier sections of this volume, subroutines CLIFT and 
DRAG arc.' also called by CTAER and AEROUT, and these various subroutines may 
be separated from one another by a number of other programs which are described 
in other volumes of this study. 

The aerodynamic input parameters are very large in number because of the 
large number of geometrically and operationally independent possibilities in 
aircraft design. In this user's manual, the .•'.nput and output quantities of 
each subroutine are listed, with their appropriate units, in alphabetic order. 
Subroutines covered in this manner are 

Table I I I, 2.1 - Subroutine AERO Input 

1 1 1. 2. 2 - Subroutine AERO Output 

111. 2. 3 - Subroutine AEROUT Input 

1 1 1. 2. 4 - Subroutine AEROUT Output 

1 1 1. 2. 5 - Surboutine CLIFT Input 

III. 2.6 - .'subroutine CLIFT Output 

1 1 1. 2. 7 - Subroutine DRAG Input 

1 1 1. 2. 8 - Subroutine DRAG Output 

1 1 1. 2. 9 - Subroutine CTAER Input 

111. 2. 10 - Subroutine CTAER Output 
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Table III. 2. 11 

111. 2. 12 

111. 2. 13 

111. 2. 14 


- Subroutine 

- Subroutine 

- Subroutine 

- Subroutine 


FIAPS Input 
FLAPS Output 
APPFLP Input 
APPFLP output 
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TABLE III. 2.1 SUBROUTINE AERO - INPUT 


VARIABLE 


DESCRIPTION 


aspect ratio 
length of tip tank, ft 
Wing span, ft 
diameter of tip tank, ft 

mean aerodynamic chord of horizontal tail, ft 

mac of vertical tail, ft 

mac of wing, ft 

average nacelle di^uneter, ft 

drag coefficient increment 

incremental effective flat plate area of fuselage, sq ft 
quarter chord sweep of wing, deg. 
incremental wetted area/wing area 
Length of tail boom, ft 
f iselage length, ft 


i.LFFC 


tusclage length minus tail boom length, ft 


' I nacelle length, ft 


uu;ar 


landing gear configuration 


■'CC NFG 


fuselage configuration 


nacelle drag estimation flag (0, 1, or 2) 
Reynolds number per foot 

relative posi .ion of horizontal tail position 


:CFAC 


factor on shift of drag divergence Mach number for super- 
critical airfoil 


fuselage wetted area 







TABLE I I I. 2.1 SUBROUTINE AERO - INPUT 
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TABLE I I I. 2. 2 SUBROUTINE AERO - OUTPUT 


VARIABLE 

DESCRIPTION 

t'BA!<l-' 

aorodynamic cleanliness coefficient 

CDNl 

skin friction coefficient of nacelle 

CFIN 

7 

skin friction coefficient at Re = 10 

CKF 

form factor of fuselage 

CKIIT 

form factor of horizontal tail 

I'KN 

form factor of nacelle 

CKTP 

form factor of tip tank 

CKVT 

form factor of vertical tail 

CKW 

form factor of wing 

CKWX 

form factor of wing 

DLTAFE 

incremental equivalent flat plate area, sq ft 

EM 

Mach numlier 

FE 

fotal equivalent flat plate area, sq ft 

FEF 

fuselage equivalent flat plate area, sq ft 

FEHT 

horizontal equivalent flat plate area, sq ft 

FF.N 

n>icelle equivalent flat plate area, sq ft 

FETF 

tip tank equivalent flat plate area, sq ft 

Fi;VT 

vertical tail equivalent flat plate area, sq ft 

FEW 

w i nq equivalent flat plate area, sq ft 

GKCD 

landing gear drag coefficient 

C.RI'E 

incremental landing gear equivalent flat plate array, 
sq ft 

SAl-4 

compressibility drag factors 

SAS, 6 

profile drag parameters 
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TABLE III. 2. 2 SUBROUTINE AERO - OUTPUT 


VARIABLE 


DESCRIPTION 


induct'd draq parameter 


SCFAX 


supercritical airfoil parameter (SCFAC) 


Oswald efficiency factor 


SWET 


total wetted area, sq ft 


SWETF 


fuselage wetted area, sq ft 


SWETH 


horizontal tail wetted area, sq ft 


SWETIC 


increment in wetted area, sq ft 


swi:tv 


vertical tail wotted area, sq ft 


SWETW 


wing wetted area, sq ft 


form factor of nacelle 
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TABLE III. 2. 3 SUBROUTINE AEROUT - INPUT 


VARIABLE 


DESCRIPTION 


array of lift coefficients 


AU'IIU) 


zero lift anqle of attack, deg 


array of Mach numbers 


aspect ratio 


CLALPH 


lift curve slope, per rad 


DLMC4 


quarter chord sweep angle, deg. 


EMCRU 


cruise Mach number 


altitude in wing spans 
= 1, calculate only low speed polar 
= 2, calculate cruise polar and low speed polar 
= 3, calculate only cruise polar 


» 4, calculate and plot botJi polars 
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TABLE III. 2.4 SUBROUTINE AEROUT - OUTPUT 
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TABLE III. 2. 5 SUBROUTINE CLIFT - INPUT 


VARIABLE 

Ai-riiA 

A1.PHLO 

AR 

CL 

CLMXLD 

CLMXTO 

DCLLD 

DCLTO 

01^04 

i;m 

HC'B 

ILIFT 


DESCRIPTION 

angle of attack, dog 
zero lift angle of attack, deg 
aspect ratio 
lift coefficient 

maximum lift coefficient in landing configuration 

maximum lift coefficient in takeoff configuration 

landing lift coefficient increment 

takeoff lift coefficient increment 

quarter chord sweep angle, deg 

free stream Mach number 

altitude measured in wing spans 

= 1, 2, 3 - compute CL from ALPHA 

= 4, 5, G - comi)ute ALPHA from CL 





i 


TABLE III. 2.6 SUBROUTINE CLIFT - OUTPUT 

VARJABLE I DESCRIPTIIoT™^ 


ALPHA aiHjlu of attack, deg. 

CL lift coefficient 

ClALI’H lift curve slope, per rad 




TABLE I I I. 2. 7 SUBROUTINE DRAG - INPUT 


VARIABLE 

DESCRIPTION 

AC DC DR 

array of ratio of win<j profile drag coefficient to 
minimum drag 

ACLS 

array of lift coefficients, dimension KWCD 

CDNAC 

nacelle drag coefficient 

CFIN 

skin friction coefficient at RE » 10 million 

CKWX 

wing form factor 

CL 

lift coefficient 

DCDAPP 

drag coefficient increment in approach configuration 

DCDLD 

drag coefficient increment in landing configuration 

DCDTO 

drag coefficient increment in takeoff configuration 

DCL.APP 

lift coefficient increment in approach configuration 

DCLLD 

lift coefficient increment in landing configuration 

IX’LTO 

lift coefficient increment in takeoff configuration 

EM 

Mach number 

GRCD 

drag coefficient of landing gear 

HOB 

altitude in wing spans 

I DRAG 

= 1 to 8, flight condition at which drag is calculated 

KNAC 

= 0, nacelle drag accounted for in engine performance 


= 1, nacelle drag part of aero drag, computed when 
engine sized 


= 2, nacelle drag part of aero drag, based on input 
nacelle size 

KVJCD 

number of points in wing profile drag table 

SAl-7 

1 to 4, compressibility drag factors 


5 and 6, profile drag parameters 


III. 2 


11 





TABLE III. 2. 7 SUBROUTINE DRAG - INPUT 


VARIABLE 

SCFAX 

SIGMAP 

SIGKLD 

SIGMTO 

VDELGA 

VDEL6L 

VDEL6T 



supercritical airfoil parameter 

empirical induced-drag correction factor, approach 
configuration 

empirical induced-drag correction factor, landing 
configuration 

empirical induced-drag correction factor, take-off 
configuration 

efficiency factor correction factor for approach flaps 
efficiency factor correction factor for landing flaps 
efficiency factor correction factor for take-off flaps 
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TABLE II I. 2. 8 SUBROUTINE DRAG - OUTPUT 








TABLE III. 2.9 SUBROUTINE CTAER - INPUT 


VARIABLE 

DESCRIPTION 

AD 

angle of attack, degrees 

ALPHLO 

zero lift angle of attack, deg 

AR 

aspect ratio 

CLALPH 

lift curve slope, per rad 

DLMC4 

quarter chord sweep angle, deg 

FPMN 

Mach number, used in fixed pitch propeller option 

H 

altitude, ft 

IDC 

reference area flag 

RWCRTX 

ratio of cruise weight to design gross weight 

SW 

wing area, sq ft 

WG 

gross weight at take-off, lb 

WGS 

wing loading, lb per ‘--q ft 
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TABLE III. 2. 10 SUBROLTINE CTAER - OUTPUT 
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TABLE III.2.U SUBROUTINE FLAPS - INPUT 


VARIABLE 

DESCRIPTION 

AAKW 

array of wing asjioct ratio 

ABFOB 

array of flap span to wing span ratios 

ABLEOB 

array of leading edge slat span to wing span ratios 

ABTEOB 

array of trailing edge device span to wing span ratios 

ACFOC 

array of flap chord to wing chord ratios 

ACROBL 

array of root chord to body length ratios 

AIX'DOT 

array of flap-induced drag coefficients 

ADCLTE 

array of flap-induced lift coefficients 

ADELCL 

array of lift coefficient increments 

ADELFD 

array of flap deflections, deg 

ADELTO 

array of optimum flap deflections for different flap typ> 

APELl 

array of flap chord correction on drag (JFLTYP )> 2) 

ADELIY 

array of flap chord correction on drag (JFLTYP £ 2) 

ADEL2 

array of flap angle correction on drag 

ADEL3 

array of flap span correction on drag 

ADELG 

array of flap angle correction of Oswald efficiency facti> 

AlAVOB 

array of wing/body diameter ratios for fuselage drag 
increment 

ALAMl 

sensitivity of maximum wing lift coefficient to aspect 
ratio 

ALA.M2 

sensitivity of maximum wing lift coefficient to thicknes:, 
to chord ratio 

ALAM3 

sensitivity of maximu-n flap lift coefficient to aspect 
ratio 

AIAM4 

ensitivity of maximum flap lift coefficient to thicknes;. 
CO chord ratio (JFLTYP ^ 2) 
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TABLE III. 2. 11 SUBROUTINE FLAPS - INPUT 


VARIABLE 

DESCRIPTION 

ALAt'l4X 

sensitivity of maximum flap lift coefficient to thick- 


ness to chord ratio (JFLTYP ^ 2) 

ALAM5 

sensitivity of maximum flap lift coefficient to flap 


chord to chord ratio (JFLTYP ^ 6) 

ALAM5X 

sensitivity of maximum flap lift coefficient to flap 


chord to chord ratio (JFLTYP = 3, 4, 5) 

ALAM5Y 

sensitivity of maximum flap lift coefficient to flap 


chord to chord ratio (JFLTYP = 1, 2) 

ALAM6 

sensitivity of maximum flap lift coefficient to flap 


angle (JFLTYP > 6) 

ALAM6X 

sensitivity of maximum flap lift coefficient to flap 


angle (JFLTYP = 3, 4, 5) 

ALAM6Y 

sensitivity of maximum flap lift coefficient to flap 


angle (JFLTYP = 1, 2) 

ALAI-17 

sensitivity of maximum flap lift coefficient to flap 


span 

ALAMIO 

sensitivity of maximum slat lift coefficient to slat 


deflection 

ALAMll 

sensitivity of maximum slat lift coefficient to slat 


span 

ALAM13 

sensitivity of maximum lift coefficient to Reynolds 


number 

ALAM14 

sensitivity of maximum lift coefficient to Mach number 

ALTFLP 

altitude at which flaps are retracted, ft 

AMNW 

array of wing Mach numbers 

AR 

wing aspect ratio 

ARDELF 

array of ratios of flap deflection to optimum deflection 


angle 

AROELL 

array of ratios of leading edge slat deflection to | 
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TADLE lit. 2. 11 SUBROUTINE FLAPS - INPUT 


VARIABLE 

DESCRIPTION 

AliNW 

iirray of Reynolds numbers for Reynolds number correction 



ASIGMA 

array of induced drag correction factors 

ASLM 

array of taper ratios in table 

ATCW 

array of thickness ratios in table 

AWCFLP 

array of weight coefficients for different flap types 

B 

wing span, ft 

BENC^B 

fraction of span blanketed by wing mounted engine 

BTEOn 

trailing edge flap span divided by wing span 

CBARW 

mean aerodynamic chord of wing, ft 

CF<.K’ 

flap chord divided by wing chord 

CUaX' 

leading edge flap chord divided by wing chord 

CKOl.W 


CKi.xn’w 

wi..g root chord, ft 


quarter chord sweep angle, deg 

KU-' 

fuselage length, ft 

II- lAP 

= 1 or 2; takeoff or approach configuration 

JFI/rYF 

= 1 to 7, plain, split, slotted, etc., flap type 

ROI»'-AX 

reference maximum lift coefficient 

RICJ'Ar 

induced drag correction factor, approach configuration 

•SIGMLD 

induced drag correction factor, landing configuration 

SIC.MTO 

induced drag correction factor, takeoff configuration 

iU>M 

wing taper ratio 

SWl' 

fuselage width, ft 


1 1 : . 2 
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TABL»: III, 2. 11 SUBROUTINE FLAPS - INPUT,' 


VARIABLE 

swi'u: 

SWI’TE 

TC 

WC.S 
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TABLE IU.2.12 SUBROUTINE FLAPS - OUTPUT 


VARIABLE 

DESCRIPTION 

CU^IXAP 

maximum lift coefficient in approach configuration 

CmXLD 

maximum lift coefficient in landing configuration 

CUIXTO 

maximum lift coefficient in takeoff configuration 

DCDAPP 

drag increment in approach configuration 

nCDLD 

drag increment in landing configuration 

DCDOTE 

drag increment due to trailing edge flap 

Dctyro 

drag increment in takeoff configuration 

DC LAPP 

lift increment in approach configuration 

DCLLD 

lift increment in landing configuration 

DCLMTE 

lift increment due to trailing edge flap 

LK’LTO 

lift increment in takeoff configuration 

Din.FD 

flap deflection, deg 

DKLTEO 

optimum flap deflection, deg 

DFI,PLD 

landing flap deflection, deg 

VDFLbA 

efficiency correction factor for approach flaps 

VDFLGL 

efficiency correction factor for landing flaps 

VDFLbT 

efficiency correction factor for takeoff flaps 

V;;tall 

stall velocity, kts 

WCFIAP 

flap weight trend coefficient 

wu:d 

weight of loading edge device, lb 

WFLJ\P 

weight of flaps, lb 

WllLDL’./ 

weight of high lift devices, lb. 
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TABLE III. 2. 13 SUBROUTINE APPFLP - INPUT (TURBOFAN) 


VARIABLE 


DESCRIPTION 


CLMXAP 


maximum lift coefficient, approach configuration 


CLMXLD 


maximum lift coefficient, landing configuration 


DELFDX 


maximum flap deflection, deg. 
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TABLE III. 2. 14 SUBROUTINE APPFLP - OUTPUT 


VARIABLE 


DESCRIPTION 


DFLAP 


approach flap deflection, deg 


= 0, computation diverges 
= 1, computation converges 
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III. 3 PROGRAMMERS MANUAL FOR AERODYNAMIC ROUTINES 

The function of the aerodynamic routines has been outlined in Volume 
III.l. Detailed flow charts ^md subroutine outlines are presented in the 
present volume. A list of inputs and outputs by the subroutine are available 
in Volume III. 2. 
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1 1 1. 3.1 Subroutine AERO - 
Flat Plate and Wetted Areas, Profile Drag 

The function of this subroutine is the computation of equivalent flat 
plate areas, wetted areas, and profile drag for the basic aircraft in the 
cruise configuration. 

Figure III. 3.1 presents a detailed flow chart of subroutine AERO. 

Calculations performed are described in Volume III.l in the following 
sections : 

III. 1.1 - Flat Plate and Wetted Areas 

1 1 1. 1.1.1- Form Factors 

111. 1.1. 2- Skin Friction Coefficient 

111. 1.1. 3- Reynold's Number Correction 

111. 1.1.4- Profile Drag Coefficient 

111. 1.1. 5- Wetted Areas and Flat Plate Areas 

11 1. 1.1. 6- Concluding Computations 

Equations (III. 1.1) to (III. 1.43) can readily be identified in the detailed 
flow chart of this subroutine 
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AERO 


o 



FIGUPE III. 3.1 - DETAILED FLOW CHART, SUBROUTINE AERO 
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AERO 


o 


ci)wi=cnN 

CI)NI=cnN 

cnri=(;i IN 

CI)VTI=CE1N 

cnnTi=i:riN 

cirnp=cnN 

RI,MC4 = I)I.M(;4*. 017453 
T= S I N ( R l>1C4 ) /COS ( R 1.MC4 ) 

YAI,i:05=n.-Sl,M)/(l.*Sl,M) 

I)l^1l’S=ATAN(T-4. /AR* ( XCPS- . 25) *YA1,E05) *RTOD 
t)l>rrCX=ATAN(T-4 . /AR* (XCTCMX- . 25) *YAU;05*RTOD 
SA2= - . 33* ( . b5-XCPS) *(!.♦. 0033* (4 . *DLMPS- 3 . *DLMTCX) ) 
RLMLE=ATAN(T^YALE05/AR) 

FK= 1 . / ( 1 . +YA1,E05/AR*4 . *SLM**2) 

SAl = ( 1 . ♦ . 0033* (4 . *Dl>lPS-3 . *DLNrrCX) ) * ( 1 . - 1 . 4*TC- . 06* (1 . -XCPS) ) -0368 
SA3= ( 1 . 5 - 2 . * FK* * 2 *S I N (RLMLE ) * * 2 ) *TC * ( . 5/ 3 . ) 

SA4= . 75*TC 
ri;f=ri:i,i*i;i,f 


Ri f=rf;li*f.lffc 


F' F' R F. = ( A LOt, 1 0 ( RF; F ) / 7 . ) * * ( - 2 . 6 ) 

RFW=RF.U*CBARW 

FWRF= ( ALOC. 1 0 ( RFW) / 7 . ) * * ( 2 . 6) 

FBMRF.= 0. 


DBOOM=SlVBM/3. 14/Hl.BM 


rf:bm=rfli*flbm 

V BMRi;= ( A LOC. 1 0 ( RE BM) / 7 . ) * * ( - 2 . 6) 


CKN= 1 . 5* (1 . ♦ . 35/ (EI.N/DBARN) ) 










AERO 


G 
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AERO 


o 


|GRCD«GRFE/SW | 



SA5=S/.5^GRCD 

GRCD=0. 




WRITE (6)CKW, 
CKVT.CIOfr.CKF, 
CKN.CKTP.CtCBM, 
DELFE 


* 

-ii 

c 

KWRITE ^ 


ll 

^ - - "■ 

-0 ^ 


^^SSsSS^Bi 


WRITE (6)CDWI 
ACLS. , I*l ,KWCD, 

c 

KWCD 



— 4 

ACDCDR. ,1=1 ,KWCD 


SAb=CKWrWRH 

SKI;Tf'=SF^SWBW 

SWhTW=:. . VSW-2. *CBARW*SWF* (1 . - . 5*SWF/B^ . 5*SLM*SWF/B)/ (1 . ♦SIM) 
SWFTV=2. *SVT 

swi;th=2.*sht 

SWr.TIODl.SWSWSW 

SW1;T=SW1;TF+SN*STIP^SWETW^SWETW^SWETH>SWETIC 

Fi:w:=sw*a)Wi*(:KW*FWRE 

FF.F--c:KF*SFn:i)FI*FFRE^DELFE^CKBM*SWBM*CDTIP*FBMRE 
Fl;VT-CKVT»SVT*C[)VTI *FVTRF 
FI;N = CKN*SN*CI)M *FNRE 
l)I.TAFF. = l)I.U:i)*SW 

FE = Fi;W*FFF*FhVT*FEHT^FEN4FETP>DLTAFE 

i;i)po=(fe-ff;w)/sw 


KNAC • |CDPO=(FE-FEW-FEN)/SW | 



RETURN 
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111. 3. 2 StibroutiM ABROUT - 
Drag Polars 


Subroutina AEROOT prints and plots low and high spaad drag polaro using 
subroutines CLIFT and DRAG to determine ai^le-of-attack or lift coefficient 
and drag coefficient. Stand^u^d system graphics routines are employed to 
prepare plots including 


Subroutine 


PLOTS 

~ 

Plot Initializer 

PLOT 

- 

Plotting Routine 

LINAXS 

- 

Linear Axis 

SYMBOL 

- 

Symbol Drawer 

SPLINE 

- 

Spline Fit 

NUMBER 

- 

Number Drawer 

RSTR 

- 

System Restoration 


A detailed flow chart for subroutine AEROUT is presented in Figure III. 3. 2. 
A discussion of this routine's function is presented in Section III. 1.2. 
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AER0UT 


o 



FIGURE III. 3. 2 “ DETAILED FLOWCHART, SUBROUTINE AEROUT 
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AER0UT 
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AER0UT 


G 


I IF(NPC .LT. 4) I 

CALL PL0TS(1O) 

CALL PL0T(O . f 




D0 130 1-1,5 
^ f 

r ’ 

D0 130 J=l,8 
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AER0UT 
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AER0OT 
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Al Hv^r ( 


CAM. I.INAXS( 1 . 7''),1 .2S.7 . 7'‘),S. .>S, . 1 2'i ,-l . 10. 1 .0, 
('I)MA.X.’, 1. 10. rDKAC) 

('AM. I.INAXS( 1 . 7S. 1 . 7 . 7S.'.. . 1 I . 10. 1 .CM, 

CI.MAX . ). lO.TLI IT) 

CAM. SYMBOL 0. 6. S. 10.0. .TUI, K 1.0. . 20) 

CALL SYMB0L(4.O.4.H').O. I .Tm,L4.0. .20) 

CALL SYMB(/»L(4.0.4. 70.0. 1 ,'l I l l.KS.O. . K') 


ni/t I .'M 1 I 


DlT* 1 l-< I- ! . 10 





AAi IH ( 1 1 \i III ( 1 

. I > 


. A '1 1 1 I 1 ••'III. 

1 1 

J 





lF(AACI.(in) .1.1. AACI,(I0-1)) 


N 


!F(ii) .i;o. 0) 


N 


1 1 1 . i 


UKIGINAL PAGE IS 
OF POOR QUALITY 































XXX. 3. 3 Subroutin* APPPLP - 
Flap implication Routine 

This subroutine controls the flap setting calculatim using stibroutine 
FLhPS to coapute lift associated with a given flap setting and XTRMHN to 
Iterate for the desired flap setting. When convergence cannot be obtained, 
flaps are set to the landing position. Refer to Section XXX. 1.3. 4 for 
additional details. A flow cheurt for subroutine APPFLP is presented in 
Figure III. 3. 3. 
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APPPLP 


4 


NPCJ-NPC 

MPC-4 

IOK-0 

DELFDX-10. 

ERROR-999. 

F-1.1 

FF-1.01 

JC-0 

JX-0 


ML 


ERRMl-ERRQR 


CALL FLAPS 


7 


c 


CLMXLD 


I 


> 


-CLMXAP 


< 


.jx.gtT^ 


I 


ABS (ERROR) .LT. 0.005. or. 


N 


I 

R 1 

I 


/ CALL ITR~ MHW(ERR0R,ERRM1>DELFDX,F .FF.JC,JXH^ 

N ^ 


< 


JC.GT.30 


DELFDX.GE.0.0 




1 DELFDX-0.0 1 

1 




L 


WRITE (6) 


FAILED TO FIND! 
FLAP SETTING 
USE LANDING 
FLAP VALUE 


CALL FLAPS 


7 


30 


NPC=NPCJ 

DFLAP=DELFDX 


-RETURN 


FIGURE III. 3. 3 

DETAILED FLOW CHART, SUBROUTINE APPFLP 
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ZZZ.3.4 SubroutiiM CZJFT • 

Wing Alo^e Lift Coefficient 

Thie subroutine conputee wing alone lift coefficient at a given angle-of- 
attack. Methodology ia described in Section ZZZ. 1.2.1 of thie report. A 
detailed flow chart of the calculations is' presented in Figure ZZZ.3.4. 
Equations ZIZ.1.54 to ZZZ. 1.65 of the referenced section can readily be 
identified in the flow chart. 
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CLIFT 


KCLGE-1.0 

DELCL-0.0 

CLMAX>2. 

RIJC4«DLMC4*.0174S 

Cl-(AR/(2.*COS(RLMC4)))**2. 

C2» 1 . - (OM*COS(RLMC4))**2. 
CLALPH-(3.1416)/(1.*SQRT(1.*C1*C2)) 


KCLGE- 1 . 005^ (.00211-. 0003* (AR-3. ) ) *C3 


J.OGr;=CLALPM* . 0174533* (ALPHA-ALPHLO) 405101 


=CLMAX 


AI,PHA=AI.I'HL04 ( (CLMAX-DELCL)/ (CLALPM* . 0174533) ) 


.L= KCl.Gi;*CI.ALPH* . 0174533* (ALPHA- ALPHLO) 4DELCL 


CL0GE=(CI,-()r.!,CI.)/KCLGE4nELCL 


>CI>1AX 


ALPHA=(CL-0ELCL)/(KCLGE*CALPH*.0174533)4ALPHL0 


? FIGURE III. 3. 4 - DETAILED FLOWCHART 

SUBROUTINE CLIFT 












IIZ.3.5 Subroutine CTABR - 
Lift and Drag in Cruise Flight 

This subroutine determines lift coefficient to support the aircraft 
in cruise flight together with the associated drag coefficient. Subroutine 
CTAER anploys 8ubr'>utine TPAJ.T to obtain atmospheric properties in the cruise 
flight condition. Subroutine CLIFT provides lift coefficient; subroutine 
AERO provides profile drag, and subroutine DRAG computes induced drag and 
ccnpressibility drag. The computations of subroutine CTAER are discussed in 
some detail in Section III. 1.3 and Equations (III. 1.71) to (III. 1.74) of 
that section can be identified in the detailed flow chart of Figiure III. 3. 5. 
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FIGUP£ I I I. 3. 5 - DETAILED FLOW CHART SUBROUTINE CTAER 












CTAER 


0 
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XXX. 3. 6 Subroutine DRAG - 


Xnduced and CompreaBiblllty Drag 

Subroutine DRAG conputes the total aircraft drag coefficient with or 
without flaps deployed. The methodology is presented In Section XXX. 1.2. 2 
of this report. Flight condition is selected through the Indicator XDRAG. 
Equations (III. 1.66) to (XXX. 1.70) can be identified In the detailed flow 
^h^u:t of Figure I IX. 3. 6. Several of the coefficients appearing in these 
equations are obtained hy interpolation in tabulated values using subroutine 
ITRLN which is the only subroutine called by DRAG. 


* 
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DRAG 



FIGURE III. 3.6 - DETAILED FLOWCHART SUBROUTINE DRAG 
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DRAG 


O 


i(EMD^((SA3^SA4*CL**2)/10.)**(l./3.)) 


WRITE (6) 


RETURN 


















DRAG 


c 


O 


I DRAG 


■S or »6 


^DELCDG»0.7^ 


C 


I DRAG 


- ^DBLCDG«0.'71 


C 


|KCDIGh=l.ll*S.S5*(HOB)-SQRT(29.8*({HOB^.02)**2)*.817)l 

L, 


KCDIGE 


I 


KwcD y~^ 



' 

|CLS=C 

i.-sigma*delcl| 

' 



n RLN ( ACLS , ACDCDR .CDCDR . KWCD) j 


CUWI=CFIN 

Cl)WI=CUWI*CDCDR 


fCD=SA5-fSA6>UELCDF>DELCDG»SA7*KCDIGE*(CL-SIGMA*DELCL)**2/VDELn 


C 


c 


ICON AC 


KNAC 

V 




- Tcd^cd^cdnac I 
- i 


[lCDNAC=o1 

RETURN 
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XII. 3.7 Subroutine FLAPS - 


Flap Characteristics Routine 

Subroutine FLAPS computes flap lift jmd drag increments for several 
flap types using the methods outlined in Section III. 1.4.1. The routine 
also has the ability to compute flap weight characteristics by the method 
of Section III. 1.4. 2. 

A detailed flow chart for subroutine FLAPS is presented in Figure ILL. 3. 7. 
The first two-and-oue*half pages of this routine are mainly concerned with 
computation of the factors VlAMl to VLAM14 (Section III. 1.4.1) which enter 
into the maximum lift coefficient computation, ihese factors are obtained 
by t£d>le look-ups using subroutine ITRLN. This is followed by computation of 
the flap efficiency factors on Pages 3 and 4 of the detailed flow chart. 

Page 5 of the flow chart primarily computes flap weights and geometrical 
characteristics. When NPC » 4, the routine transfers directly to the bottom 
of Page 6 at Statement 150 where allowable approach flap characteristics 
are determined. 
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FLAPS 


o 



RETURN - 


Rl.MC4»DLMC4*. 017453 


* NPC 


= 2 =3 


DCLMTE 


iDCLhfTE-ADCLTE ( JFLTYP] 


OELTEO 


IDELTFO KLTO (JFLTYP) 


DC DOTE 


pc DOTE=ADCDOT (JFLTYP) 


WCFLAP 


WCFLAP=AWCFLP (JFLTYP] 


CALL ITRLN(/URW.ALAM1,AR . V LAM 1 , N ARW ) 


[call I'l rln (A Tt:w , al/\m: ,rc , vi,am2 , nicw) 


[CALL I'RLN (AAKW , AL/\M5 , AK, V LAM5 , NARW j 


JFLTYP 


[CALL 1TRLN(ATCW,ALAM4X,TC,VLAM4,NTCW) 


JFLTYP 


CALL ITRLN(ATCW,ALAM4,TC,VLAM4,NTCW) 


JFT.TYP 


CALL 1 rKl,N(ACFOC.ALAM5Y,Cl()(:,VI,AM;. .NCIOC)| 


JFLTYr 


■ f 1. . = ~^ fcAUnmn7r?ACF0^ 


JFLTYP 


CALL n'RLN(ACF0C,AI.AM5.CF'(K:,VLAMr.,\CI (X;) 


JFLTYP 


CAI.L ITRLN(ADi;U-D,ALAM(>Y.Di;i,l I),VL\M(.,N1)1,1,1'D) 
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FIGURE III. 3. 7 

detailed flowchart subroutine Fi.AI'S 




FLAPS 


and ■$ 


O 


CALL ITRLN ( ADELFD , ALAM6X , DELFD , VLAM6 , NDELFD) 
















FU^S 


o 



* MAxifflum Lift 
Coefficient 


Ql-WGS/CLMX 
P0»P0*144. 
SMN*SQRT(Ql/0.7/P0) 
VK«SMN*SA/ 1.689 


c 


VK 


> 


<0 


SMN=VK*1.689/SA 

Ql=0.7*P0*SMN*SMN 

CLMX=WGS/Q1 

RNW=CBARW*VK* 1 . 689/XKV 
RNW=RNW/ 100000. 


I 


Hvk«o. 


fcATj^TmUAMNW^^ 


|CLM\X=(R(:i>lAX*V'LAMl*VLAM2^DCl>rrH*VLAM:<*VLAM4*VLAM5*VLAM6*VLAM7* 

VL,\M8>I)(:LMI.I;*VLAM9*VLAM10*VI,AM1PVLAM12)*VLAM13*VIJU^I4*DELCLF| 

|i;rror=(:i,max-c.,mx 


KWRITl-: 




WRrn;(6)HKR0R,LRRMl, 

VK,CLMAX,CI.MX,VLAM13, 

VLAM14,RNW,SMN,JC 


c 


ABS (ERROR) 


2) 0 


I 




c 


.IC 


>10 


.n-LTYP ^ 


= 2 


CALL lTRLN(ACFOC,ADELlY,CFOC,VnnLl .NCFOC) 


C 


JFLTYP 


— -»•> — fTAu" 


ITRLN(ACF0C,A1)LL1 .Cl'OC . VI)I;I,I .NCIOC) 


I Rl)FiLF=l)Ll.FD/nF.LT(-:0 1 


^ CAM, ITRI,N (ARI)|;LF 

i ._ , 

I CALI, BlV(VI)lil.3,Hri:UB,SI>1,ABFOB.ASLM,AUl:L3.NBI'OB,NSl>l,NliR)| 


c 


NIiR 
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-0 



FLAPS 


O 


TSWPFH-TAN (SWPTE) ♦4. *CFOC* (TAN (RLMC4) -TAN (SWPTE) )/3. 


SWPFHL»ATAN(TSWPFH 

NERPT*19 

VDEL4=C0S(SWPFHL) 

NERPT=20 

VDEL5«1.0-FW0B 


CALL I I K LN ( AOf • LFU . ADE L6 , DE LFD , VDE L6 , NDELFD) 


Cro=SAA5>SA6*CFIN 

CLMNHL=RCLMAX*VLAM1 *VLAM2*VLAM1 3*VLAM14 
CDACLM=CDO^DCDOTE*VDELl*VDEL2*VDEL3*VDEL4*VDELS 
♦(CLMHL**2)/(VDEL6/SA7) 

DCL=DCLFfrE*VLAM3*VLAM4*VLAM5*VLAM6*VLAM7*VLAM8*VLAM13*V!JVM14 

DCD=DCD0TE*VDEL1*VDEL2*VDEL5*VDEL4*VDEL5 


CLMXTO=CLMAX 

DCLTO=DCL 

DCDTO-DCD 

VDEL6T=VDEL6 


CALL rrPLN(ADELFD,ASIGMA,DELFD,SIGMT0,15) 


WRITE (6)CLMXT0, 

VK.DELFD.DELLED, 

DCLTO.DCDTO 




WRITE (6)CLMXT0, 

VK.DELFD.DELLED, 

DCLTO.DCDTO 


CLMXD=CIWAX 

D!'U’l.l)=f)i;LFD 

DCLLD=I)CL 

DCDLD^DCD 

Vl)|-.),Ol,=Vl)EL6 


CAM, nUI.N(AI)l.LFD,ASIGMA.OELFD,SIGMLD,15) 
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LAPS 


WRITE(6)CLMXMD, 
VK.DELn.DF.LLED, 
DCLLD.D'jDL 




.-SLM)*FWOB- 
( 1 . -SLM) * BTEOB) * (BTEOB- FWOB- BENGOB) ) / . 2 


SFLAP»0. 


SA I L» { (CFOC*CRCLW* B) * (SLM^ 1 . - ( 1 . -SLM) * BTEOB) * ( 1 . -BTEOB) ) / 2 . 
SLE= ( (CLEOC*CRCLW*B) * (1 . - (DBALE/B) * (SLM- 1 . ) ♦SLM) *BLEOB) /2 . 
SBOX* ( 1 . O-CFOC-CLEOC) *CRCLW*B* { 1 . O^SLM) /2 . 

SENG= (CFOC*CRCLW*B* (1 . 25-FWOB^SLM* (FWOB^O. 75) j *BENGOB) /2. 
STE=SFLAP/3. ♦SENG^O. 15*SAIL 




RMO=PO/(53.32*G*TO) 

VSTALL*.59 

VFLAP=1.8*VSTALL 


WFLAP»WCFLAP* (VFLAP/ 100. ) ** 2 . 
*SFLAP*FLAPN**(-.5) 




*SFLAP*FLAPN*.5 


WFLAP=WCFLAP* (VFLAP/ 100. **2. 38 
*SFLAP**1. 19/(FLAPN**.595] 


WF’ LAP= WC F LAP * SF LAP * (V I- LAP **2.195) 
/45180. 


W I- LAP = WC F LAP * S F LAP * 0 . 369 
*VFLAP**0.2733 


WLED-3.28*SL!;**1.13 

WHLDIiV=WFLAP>WLED 


WRm:(6) 
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FLAPS 


G 













FLAPf. 


■3^ and »S 


O 


CALL ITRLN (ADELFD . A LAM6X , DE LFO , VLAM6 , NDELFD 


CALL ITRLN (ADELFD, ALAM6,DELFD,VLAM6, NDELFD 




CLMAX=(RC 

VLAM7*VLAM8>DCLMLE*VLAM9*VLAM10*VLAM11*VLAM12) 

*VLAM13*VLAM14+DELCLF 

RDELF=DELFD/DEMEO 


CALL ITRLN (ARDELF , ADEL2 , RDELF , VDEL2 .NRDELF 


CALL ITRLN (ADELFD, ADE L6 , DE LFD , VDE L6 , NDELFD) 


DCL=DCLMTE*VLAM3*VLAM4*VLAM5*VLAM6*VLAM7*VLAM8*VLAM13*VLAM14 

DCD=DCD0TE*VDEL1*VDEL2*VDEL3*VDEL4*VDEL5 

CLMXAP=CLMAX 

DCLAPP=DCL 

DCDAPP=DCC 

VEL6A=VDEL6 


ALL ITRLN (ADELFD, ASIGMA,DELFD,SIGMAP, 15) 
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